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On  the  theoretical  side,  a  new  highly  scalable  code  for  multi-configurational  self- 
consistent  field  (MCSCF)  calculations  has  been  developed.  This  code  is  in  the  latest 
GAMESS  release.  A  parallel  code  for  second  order  unrestricted  open  shelf  perturbation 
theory  (UMP2)  gradients  has  been  developed.  The  derivation  for  the  Z-averaged  open 
shell  second  order  perturbation  theory  gradients  has  been  completed,  and  a  paper 
describing  this  derivation  has  been  published.  The  parallel  code  for  this  method  has  now 
been  completed.  A  new  second  order  multi-reference  perturbation  theory  (based  on  the 
generalized  Van  Vleck  perturbation  theory),  GVVPT2  has  been  implemented  in 
GAMESS.  A  preliminary  code  for  the  GWPT2  gradients  has  been  written,  and  a  more 
robust  code  is  in  progress.  The  MCSCF  analytic  Hessian  code  has  now  been  completed. 
Initial  scaling  tests  suggest  that  this  code,  now  available  in  GAMESS,  will  scale  very 
well  with  the  number  of  processors.  A  review  paper  on  second  order  perturbation  theory 
gradients  has  been  published.  A  new,  very  effective  coarse-grained  parallel  method, 
called  GDDI,  has  been  developed  and  implemented.  This  method  has  already  proved  to 
be  very  effective  for  numerical  gradient  and  Hessian  calculations,  since  each  energy  or 
energy  +  gradient  calculation  is  completely  independent  of  each  other.  Other 
developments  include  the  implementation  of  a  new  density  functional  theory  (DFT)- 
based  effective  fragment  potential  (EFP).  method,  the  development  of  molecular 
dynamics  and  Monte  Carlo  methods  to  facilitate  the  study  of  solvation  and  liquid 
behavior,  and  the  development  and  implementation  of  several  MCSCF.-based  approaches 
to  spin-orbit  coupling.  With  regard  to  applications,  considerable  progress' has  been  made 
in  our  understanding  of  the  mechanisms  for  formation  .of  POSS  (polyhedral  oligomeric  ; ,  > 
silsesquioxanes)  and  their  titanium  analogs.  A  paper  on  ionic  liquids  based  on  ;  •' . 
triazolium  cation  has  been  published.  A  second  paper  on  tetrazolium  has  been  submitted, 
and  a  paper  on  pentazolium  cation  is  in  draft  form.  All  of  these  calculations,  have  been 
performed  at  a  reasonable  level  of  electronic  structure  theory  (MP2  and  G2  for  heats  of  # 
formation).  Studies  of  multiple  ion  pairs  have  been  initiated  using  our  effective  fragment 
potential  (EFP),  followed  by  re-optimizations  with  MP2.  ;  .  ..  .  ,  .  ;  . 


4.  Accomplishments/New  Findings:  Describe  research  highlights,  their 
significance  to  the  field,  their  relationship  to  the  original  goals,  their  relevance  to 
the  AF’s  mission,  and  their  potential  applications  to  AF.  and  diviliant technology 
challenges. 

A.  THEORY/MODEL  DEVELOPMENT  :  • r  ' 


Parallel  Developments.  A  highly  scalable  code  for  frozen  core  second  order  perturbation 
theory  gradients  for  closed  shell  molecules  has  been  developed  and  is  on-line  for  general 
use  at  several  MSRC  systems.  This  new  code  uses  the  Distributed  Data  Interface  (DDI) 


so  that  the  large  arrays  do  not  have  to  be  replicated.  This  means  that  much  larger 
molecules  can  be  studied  using  geometries  based  on  correlated  wavefunctions.  The  code 
scales  very  well  for  as  many  as  512  nodes  and  has  already  been  applied  to  several 
challenging  compounds.  The  development  of  analogous  codes  for  molecules  with 
unpaired  electrons,  using  both  restricted  and  unrestricted  wavefunctions  has  been 
completed.  A  general  overview  paper  on  second  order  perturbation  theory  derivatives  has 
been  published.  A  paper  on  a  DDI  implementation  of  the  SCF  method  has  also  been 
published.  Improved  parallel  MCSCF  codes  have  been  completed,  and  a  distributed 
parallel  full  Cl  code  has  been  implemented  and  a  paper  published.  A  scalable  MCSCF 
Hessian  code  has  been  developed  and  implemented.  The  multi-reference  perturbation 
theory  code  (MCQDPT)  has  now  been  implemented  under  DDI.  The  effective  fragment 
potential  (EFP)  code  has  been  made  scalable  under  DDI.  In  an  important  related 
development,  parallel  GAMESS  now  runs  on  clusters  of  PCs  and  Macs  (running  Linux) 
and  high  end  workstations.  This  is  more  challenging  than  self-contained  massively 
parallel  computers,  since  the  overhead  due  to  inter-node  Communications  is  more 
complex.  This  issue  has  been  solved  to  some  extent  by  using  a  Gigabit  Ethernet  switch 
with  large  data  packets.  We  are  also  exploring  alternative  communications  solutions, 
such  as  Myrinet,  SCI,  and  Infiniband.  These  developments  and  their  applications  have  ■  : 
been  enhanced  by  the  construction  of  a  new  32-node,  128-CPU  IBM  Power3+  cluster 
using  funds  provided  by  a  DURIP  grant.  As  for  all  other  GAMESS  developments,  we 
will  make  our  experiences  in  developing  scalable  clusters  available  to  all  users. 

Condensed  Phase  Methods.  We  have  already  shown  that  our  effective  fragment  potential 
(EFP)  method  for  solvation  is  excellent  for  water,  in  a  variety  of  applications.  We  are 
now  working  on  extending  the  capabilities  of  the  method  in  several  ways.  We  are 
exploring  several  alternative  approaches  for  incorporating  dispersion  and  other  higher 
order  terms  into  the  method.  Such  terms  are  particularly  important  for  nonpolar  solvents. 
We  are  also  in  the  process  of  extending  the'iribdel  so  that  it  is- equally  applifcable  and,  in¬ 

accurate  for  any  solvent.  Key  to  the’  success  is  the  derivation  and  implementation  of 
general  expressions  for  charge  transfer  and  dispersion  contributions  -that  contain  no  fitted 
parameters.  Now  that  these  new  terms  have  been  implemented,  an  EFP  can  be  generated 
for  essentially  any  species,  in  particular  ionic  liquids,__A.paper  on  the  dispersion  term  has 
been  published,  and  a  paper  on  the  charge  transfer  termhas  just  been  accepted.  The 
derivation  of  an  expression  for  the  analytic  gradient  for  the  EFT-ab  initio  interaction  term 
is  in  progress.  A  very  important  new  development  (by  co-worker  Jan  Jensen)  is  a  new 
method  for  using  EFPs  across  covalent  bonds.  This  will  facilitate  the  representation  of 
large  substituents,  as  well  as  the  treatment  of  large  biological  molecules.  We  have  also 
interfaced  the  EFP  method  with  two  continuum  methods,  the  simple  Onsager  reaction 
field  and  the  more  sophisticated  polarizable  continuum  -model  (PCM). 

As  the  number  of  solvent  molecules  in  the  system  increases,  the  number  of  configurations 
to  be  considered  increases  rapidly,  and  traditional  small  molecule  geometry  optimization 
methods  are  not  effective.  We  are  therefore  developing  both  molecular  dynamics  and 
Monte  Carlo  simulation  codes  so  that  the  configurational  space  can  be  probed  more 
effectively,  not  only  for  minima,  but  for  transition  states  and  reaction  paths  as  well.  This 
is  also  an  important  first  step  in  the  development  of  methods  based  oh  the  EFP  model  for 


this  is  a  very  efficient  approach,  its  accuracy  is  limited.  However,  there  is  an 
intermediate  method,  in  which  only  those  matrix  elements  connecting  orbitals  in  the 
active  space  with  each  other  are  neglected.  We  have  shown  that  in  this  case  excellent 
agreement  with  the  all  electron  method  is  obtained,  without  any  need  for  fitted 
parameters.  We  have  also  implemented  the  RESC  (relativistic  elimination  of  the  small 
component)  method  for  calculating  spin-orbit  coupling  interactions,  complete  with 
analytic  gradients. 

The  other  source  of  non-adiabatic  interactions  is  derivative  (vibronic  coupling).  We  have 
now  derived  and  implemented  analytic  CASSCF  Hessians  (energy  second  derivatives), 
and  the  analogous  derivative  coupling  terms  have  been  derived  and  coded.  This 
functionality  will  be  essential  for  studying  excited  states  where  surface  crossings  are 
common. 


Density  functional  theory  (DFT)  has  become  extraordinarily  popular  in  recent  years, 
because  it  has  the  same  computational  requirements  as  Hartree-Fock  theory,  but  often  has 
a  level  of  accuracy  that  is  comparable  to  second  order  perturbation  theory.  We  have 
developed  and  implemented  a  gridless  approach  to  DFT,  by  using  the  resolution  of  the 
identity  to  evaluate  the  complicated  integrals.  This  procedure  requires  the  use  of 
auxiliary  basis  sets  to  resolve  the  identity.  After  extensive  analysis  of  the  convergence  of 
several  molecular  properties  as  a  function  of  basis  set,  we  have  recommended  a  series  of 
auxiliary  basis  sets  for  use  with  the  first  row  in  the  periodic  table.  The  next  step  will  be 
to  extend  this  basis  set  development  to  heavier  elements.  We  have  also  implemented  a 
more  traditional  grid-based  DFT  method,  including  all  of  the  most  popular  functionals.^  vi 


Another  approach  that  is  very  important  for  large  molecules,  especially  molecules  that 
contain  heavier  elements,  is  the  effective  core  potential  (ECP)  method.  This  method 
replaces  inner  shell  electrons  by  model  potentials,  making  the.calculafions  much  more 
efficient.  These  potentials  include  relativistic  corrections  for  the  heavier  elements.  To 
make  these  calculations  more  useful,  we  have '.developed much  more efficient 'methods..;  * 
for  obtaining  the  first  and  second  energy  derivatives;  that. is,,  the.  gradients  and  hessians. 
The  increase  in  speed  is  a  factor  of  5-10,  and  the  method  has  been  implemented  to  take 
advantage  of  parallel  computers.  Very  recently,  our  collaborator  Mariusz  Klobukowski  • 
has  implemented  model  core  potentials  (MGP)  into  GAMESS.  Because  MCPs  have  the  . 
correct  nodal  characteristics,  they  are  more  accurate  for  many  properties,  such  as  spin- 
orbit  coupling  matrix  elements  and  relativistic  effects  for  heavy. elements.  .  ;  .  ;  .  ^  ,mv 


As  calculations  on  larger  and  larger  species  are  made  feasible  by  advances  in  parallel 
programs  and  other  innovative  algorithms,  it  becomes  increasingly  important  to  develop 
graphical  methods  for  interpreting  the  results  of  complex  calculations.  This  has  led  to,  the 
continuing  development  of  MacMolPlt,  a  graphical  interface  for  GAMESS.  MacMolPlt 
provides  the  ability  to  visualize  (and  rotate  in  real  time)  orbitals,  wavefunctions,  electron 
densities  and  electrostatic  potentials  in  two  and  three  dimensions'.  We  ate  now  irfthe 
process  of  developing  a  GUI  (graphical  user  interface)  that  will  make  the  construction  of 
large  and  complex  species  much  easier. 


We  have  also  been  exploring  alternatives  and  extensions  to  the  popular  G2  and  G3 
methods  developed  by  Pople  and  co-workers  for  the  accurate  prediction  of  such 
theormodynamic  properties  as  heats  of  formation,  ionization  potentials  and  electron 
affinities.  One  limitation  of  these  methods  is  that  they  are  applicable  only  to  species  that 
are  adequately  described  by  single  configuration  wavefunctions.  This  eliminates  many 
transition  metal  compounds,  diradicals  and  most  transition  states.  We  have  therefore 
initiated,  in  collaboration  with  the  Radom  group,  the  development  of  multi-reference 
analogs  of  the  G2  and  G3  methods.  The  multi-reference  methods  are  based  on  CASSCF 
wavefunctions,  instead  of  Hartree-Fock,  followed  by  multiTreference  perturbation  theory 
(instead  of  MP2)  and  finally  multi -reference  Cl.  The  initial  paper  describing  this  method 
has  been  published. 


B.  APPLICATIONS 

POSS  Compounds.  Polyhedral  oligomeric  silsesquioxanes  (POSS)  have  a  wide  Variety  of 
important  applications,  including  lubricants  and  coatings.  They- are. also  resistant  to 
extreme  environments.  Consequently,  there  is  considerable  interest  in  these  species  in 
industry,  universities  and  Air  Force  Laboratories.  We  have  therefore  initiated  a  detailed 
study  of  POSS  compounds,  with  primary  emphasis  on  the  mechanisms  by  which  these 
species  form,  in  both  the  gas  phase,  and  in  solution  and  the  presence  of  catalysts.  In  the 
first  series  of  calculations,  we  examined  the  three  stepwise  hydrolysis  steps  of 
trichlorosilane,  followed  by  the  condensation  of  the  product  trihydroxysilane.  All  four 
reactions  are  predicted  to  have  large  barrier  heights  in  the  gas  phase^  but  the  addition  of 
just  one  extra  water  molecule  is  sufficient  to  reduce  all  barriers  to  zero  or  nearly  zero, 
except  for  the  first  hydrolysis  step.  In  the  latter,  there  is  still  a  residual  barrier  of  almost 
10  kcal/mol.  So,  the  next  step  will  be  to  examine  the  effect  of  additional  water 
molecules,  especially  on  this  first  step..  «•-.  -:v ;  surface,  it  was  discovered  that  there.  are  v.- 

Calculations  have  also  been  completed  on  the  next  steps  in  the  mechanism,  in  which  the 
initial  condensation  products,  disiloxanes,  are  further  condensed  to  the;  ring  compounds^ 
D3  and  D4.  These  results  are  consistent  with  those  reported  above,  in'  that  initially  high'  ’ 
barriers  are  reduced  to  nearly  zero  by  an  additional  water  molecule.  The  most  recent 
effort  has  been  the  study  of  the  addition  of  more  water  molecules,  which  lower  barriers 
even  further,  and  analysis  of  substituent  effects  on  the  reaction  mechanism  and 
energetics.  Most  recently,  we  have  comp  leted-the.  first  systematic  studyofthe .  competing  . 
mechanisms  for  the  entire  formation  of  a  POSS,  Tg,  from  the  starting  siloxanes.  This 
work  has  now  been  published.  : t-r.—hedra!  TAf/  A.,.  :...  L. .. . 

Although  our  main  focus  is  on  Si-based  POSS;  the.  Ti  analogs  are  also  of  interest;  We/, 
have  therefore  initiated  a  series  of  calculations  on  the  Ti-POSS  compounds  to  study  their 
properties.  The  first  paper  on  mixed  Si-Ti  POSS  has  been  published.  These  mixed  Si-Ti 
POSS  also  have  potential  as  catalysts,  and  a  paper  of  the  efficacy  Of  these  potential 
catalysts  has  been  published. 


Three  dimensional  cage  compounds,  including  zeolites,  have  been  of  interest  as  a 
possible  means  for  separating  small  gas  molecules.  In  collaboration  with  Dr.  Shawn 
Phillips  (AFRL-Edwards)  we  have  therefore  completed  a  series  of  calculations  in  which 
the  potential  energy  surfaces  for  passing  N2  and  O2  through  the  faces  of  Tn,  n  =  8, 10,  12. 
For  Tg,  both  the  energy  barriers  and  the  endothermicities  are  very  large  for  both  N2  and 
O2.  In  both  cases,  electron  correlation  is  essential  for  an  adequate  estimate  of  the 
energetics.  For  N2,  the  barrier  height  is  on  the  same  order  as  the  SiO  bond  strength,  while 
the  O2  barrier  is  much  smaller.  Despite  these  observations,  the  SiO  bond  does  not  break 
upon  entry  of  either  molecule,  and  the  species  with  the  gas  molecules  are  minima  on  their 
respective  potential  energy  surfaces.  As  the  size  of  the  cages  increases,  the  barriers  and 
endothermicities  decrease  as  one  would  expect.  A  paper  on  this  study  has  been 
published.  We  are  now  exploring  similar  possibilities  for  storing  H2  in  POSS  cages. 
Preliminary  calculations  on  Ts  and  its  Ti  and  Zr  analogs  suggest  that  it  is  easier  to  insert 
H2  into  these  cages,  but  there  are  still  substantial  energy  barriers.. 


Fuels  with  High  N-Content.  There  has  been  considerable  interest  in  the  last  several  years  v  '  ■  - . 

in  the  potential  for  cubic  molecules  as  high  energy  fuels.  Cubane  itself  has  been 

considered,  and  so  far  rejected  because  of  the  complexity  of  its  synthesis.  The  octksila  ■  ’ T ’  . 

analog  of  cubane  has  been  prepared.  Unlike  the  carbon  system,  the  cubic  structure  is  the 

most  stable  SigHg  isomer.  This  makes  it  potentially  very  interesting,  since  it  is  still  rather 

high  in  energy.  However,  the  heavy  mass  of  Si  precludes  any  viability  of  such  species  as 

high  energy  fuels.  , 

On  the  other  hand,  cubic  Ng  is  a  light,  very  high  energy  fuel  that  has  been  the  subject  of 

considerable  attention.  Extensive  calculations,  however,  illustrate  that  the  upper  limit  of  . 

the  barrier  separating  this  species  from  four  nitrogen  molecules  is  20  kcal/mol.  It  is, 

therefore,  not  a  viable  high  energy  fuel..  In  the  process  of  studying  the  cubic  Ng  pbtferitial  * 

energy  surface,  it  was  discovered  that  there  are  three iN^isomers  that  hr&muchdoY^inUev*  h>. 

energy  than  the  cube,  yet  still  much  higher  in  energy  than  molecular  nitrogen,  The,  : . 

potential  energy. surfaces  for  these  isomersi have. now -beenstudiedlm detail J?TwbCof, these'.! !•.>•... 

have  barriers  on  the  order  of  15-20  kcal/mol(too  small),  but  the  third  has  a  barrier  of 

nearly  30  kcal/mol.  So,  this  isomer  may  be  a  viable  synthetic  target. 


Related  to  the  foregoing,  much  attention  has  been  paid  recently  to  N4  as  a  potential  fuel. 

Many  references  are  made  to  a  barrier  on  the  order  of  50-65  kcal/mol,  even  though  it  has 

been  shown  by  Yarkony  that  non-adiabatic  interactions  reduce  this  barrier  to  less  than  4Q, _ 

kcal/mol.  For  the  first  time,  we  have  mapped  out  the  potential  energy  surface  that  leads 
from  tetrahedral  N4,  using  second  order  perturbation  theory.  One  also  must  consider  the ; ,  =. 
possibility  that  two  N4  molecules  could  interact  and  form  molecular  nitrogen.  This 
possibility  is  being  explored  in  this  laboratory. 

Rob  Schmitt  at  SRI  proposed  an  intriguing  series  of  potential  high  energy  compounds  ; 
that  we  have  been  exploring  computationally.  The  procedure  is  to  use  a  modified  version 
of  G2  to  predict  the  heats  of  formation  and  then  to  use  these  predictions  to  assess  the 
specific  impulse.  For  the  most  promising  species,  further  studies  are  performed  in  which 


the  stability  of  these  compounds  with  respect  to  unimolecular  decomposition, 
bimolecular  interactions  and  attack  by  environmental  species  such  as  water. 

Christe  and  co-workers  at  AFRL  (Edwards)  have  isolated  the  first  new  all-nitrogen 
species  in  nearly  a  century.  With  the  counter  ion,  this  species  is  [AsF6"][N5+],  We  have 
now  studied  the  potential  energy  surface  of  this  species  as  it  forms  AsF5  +  FN5,  and 
subsequently  the  pathway  leading  from  FN5  to  FN3  +  N2.  A  paper  describing  this  work 
has  recently  been  published.  Another  paper  on  the  possibility  of  making  crystals  of 
N5VN5",  has  been  published.  Both  of  these  latter  papers  have  been  done  in  collaboration 
with  colleagues  at  AFRL- Edwards.  . 

Cryogenics.  There  has  been  considerable  experimental  and  theoretical  interest  in  metal- 
doped  solid  hydrogen.  We  have  an  ongoing  series  of  investigations  of  the  low-lying  •; 
electronic  states  of  BH2  system,  to  determine  the  energetics  by  which  the  metal  is.  held  as 
a  weakly  bound  species  in  the  matrix  and  the  ease  with  which  two  metal  atoms  find -each 
other.  Similar  studies  are  under  way  on  the  B2-H2  system  and  the  A1  analogs. 

The  desired  reaction  product  for  burning  Al-doped  H2  is  AI2O3,  since  it  is  the 
thermodynamically  most  stable  aluminum  oxide.  However,  it  is  unknown  how  or  when 
this  species  forms  in  the  process.  We  have  therefore  initiated  an  extensive  study  of  the 
potential  energy  surface  that  leads  from  A1  and  O2  to  various  oxides.  A  paper  on  the 
reaction  of  A1  with  O2  to  from  A102,  and  the  subsequent  formation  of  A1203  has  now 
been  published. 

Ionic  Liquids.  A  major  effort  in  the  study  of  ionic  liquids  has  been  initiated. 
Computations  on  the  electronic  structures  of  substituted  triazplium,  tetrazqlium  and 
pentazolium  cations  are  ongoing,  in  an  attempt  to  assess  which  substituents  best  stabilize 
these  cations  and  enhance  desirable  properties;.; ,  HeatS/Offormationqf  these  substituted 
cations  have  been  predicted  using  G2  theory, and  i^odesmic  reactions.  Calculations  on 
each  of  these  cations  paired  with  NO2',  CKV,  and  dinitramide  have  also  been  performed. 
These  calculations  suggest  than  a  single  pair  of  ions  frequently  undergoes  a  barrier-less 
proton  transfer  to  form  a  neutral  pair,  although  a  small  number  of  stable  ion  pairs  has  ~  ~ 
been  found.  This  suggests  that  one  needs  to  increase  the  number  of  dimer  pairs, 
approaching  a  unit  cell  and  ultimately  the  bulk.  This  can  be  accomplished  using  the  new 
EFP2  method,  now  that  dispersion  and  charge  transfer  terms  have  been  derived  and 
implemented.  We  are  now  in  the  process  of  using  the  EFP2  method  to  predict  structures 
and  energetics  for  multiple  ion  pairs.  For  small  systems,  the  geometries  will  be  re¬ 
optimized  using  MP2,  so  that  the  EFP2  method  can  be  carefully  calibrated.  For  larger 
systems  and  for  simulating  bulk  properties,  the  EFP2  method  will  be  used  on  its  own. 
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